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Summary. — The scale invariant Petrov classification of the Weyl tensor is linked 
to the scale invariant combination of the Kasner index constraints, and the Lifshitz- 
Khalatnikov Kasner index parametrization scheme turns out to be a natural way of 
adapting to this symmetry, while hiding the permutation symmetry that is instead 
made manifest by the Misner parametrization scheme. While not so interesting for 
the Kasner spacetime by itself, it gives a geometrical meaning to the famous Kasner 
map transitioning between Kasner epochs and Kasner eras, equivalently bouncing 
between curvature walls, in the BLK-Mixmaster dynamics exhibited by spatially 
homogeneous cosmologies approaching the initial cosmological singularity and the 
inhomogeneous generalization of this dynamics. 
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1. Introduction 

It is amazing that after years of people using and working with a piece of well known 
mathematics, there is still the possibility of realizing something new about it that some- 
how escaped notice for so long. While not earth shaking, the present remarks about the 
connections between the Weyl tensor and the Kasner index parametrization are certainly 
interesting. Discovered serendipitously by playing around with the speciality index [1], 
a little thought uncovered some cute mathematical structure [2] , aspects of which have 
played a key role in the lore of general relativity in the past half century. 

The speciality index S* is a dimensionless curvature invariant which indicates the 
the special or non-special Petrov algebraic behavior of the Weyl curvature tensor of 
a spacetime which is not too degenerate in the Petrov classification of such tensors. 
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This quantity, originally introduced and used in numerical relativity in the study of 
a black hole radiating gravitational waves [1,3], also found application to the Petrov 
classification of perturbed spacctimcs [4.5], to wave extraction in numerical relativity [6] 
and in cosmology [7]. The Kasner spacetime itself has been used as an example for 
discussing the determination of quasi-Kinnersley frames [8]. 

The l-paramcter family of Kasner spacetimes [9] is most simply represented in terms 
of the Kasner indices, three exponents satisfying two identities, and which describe the 
scalefree part of the eigenvalues of the extrinsic curvature tensor for the spatial geome- 
try. The spatial metric isometry classes correspond to the equivalence classes of Kasner 
triplets under permutations of the three Kasner axes, leading to a natural action of 
the symmetric group S3 of permutations of the integers (1,2,3). The Petrov classes 
characterizing the Wcyl curvature tensor involve the scalefree part of the eigenvalues 
of a complex matrix arising from the self-dual combination of the Weyl tensor and its 
dual. These classes correspond to equivalence classes modulo the same group action 
of S^. It is not surprising that there should be a nice relationship between the two 
mathematical structures, and it is this relationship which gives added meaning to the 
Lifshitz-Khalatnikov parametrization of the Kasner indices, introduced in order to de- 
scribe collisions with ciirvature potential walls in asymptotic gravitational dynamics in a 
simple way, now seen in terms of simple transitions of the Weyl tensor eigenvalue ratios. 

To summarize the key points which the following discussion illuminates in detail: 

1. As described in Ref. 10, the Weyl tensor Ca^-ys {ct,/3 = 0, 1,2, 3) can be repack- 
aged as a complex self-dual tensor Cai3-,s + i*Coif)'y5 which in turn with a choice 
of unit timelike vector field u" can be repackaged as a tracefree symmetric spatial 
tensor Q^b = —{Eab + iBab) {a,b= 1, 2, 3), which is a complex combination of the 
electric and magnetic parts of the Weyl tensor with respect to that vector field, 
and whose eigenvalues Aq therefore sum to zero: Ai -|- A2 -|- A3 = 0, leaving only 
two independent eigenvalues, and modulo rescaling, only one independent complex 
number to describe its scalefree part (at a given spacetime point). Any one of the 
6 possible ratios Hab = ^a/\ of the three eigenvalues can be taken to parametrize 
this degree of freedom, leading to the action of the symmetric group Sz on the 
resulting parametrizations. The classification of the eigenvalues (modulo rescaling) 
is the Petrov classification. 

2. The family of Bianchi type I Petrov type I (and D) vacuum Kasner spacetimes is 
parametrized by three Kasner indices Pa which are the eigenvalues of the scalefree 
part of the extrinsic curvature tensor K°'i, satisfying two identities whose scalefree 
combination is simply 1 / pi + 1 / p2+^ / Pd, = 0, and hence their real ratios Uab = Pa/Pb 
(six possibilities, with ^3 again acting on the choices) parametrize the scalefree part 
of these eigenvalues, only one of which is independent due to the second independent 
identity. 

3. For Kasner, the Weyl tensor ratios \a/\ agree with the Kasner index ratios 
{l/pa)/{^/pb) — Pb/Pa'- IJ'ab — Uba, thus tying the Kasner map at the heart of 
BKL-Mixmaster dynamics directly to simple transitions in the scale invariant part 
of the Weyl tensor [2] . 

After expanding on these points, we review the role of the Lifshitz-Khalatnikov 
parametrizations Uab of the Kasner indices in the BLK-Mixmaster dynamics [11-18,20- 
22]. 
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2. The Weyl tensor 

Consider a timelike congruence of world lines with unit tangent vector u" satisfying 
Uau" = — 1. u" can be interpreted as the 4- velocity field of an observer family and 
used to measure tensor fields in the 1 + 3 spacetime splitting which follows from the 
orthogonal decomposition of each tangent space into a local time direction and a local 
rest space. This description of spacetime geometry is often referred to as gravitoelectro- 
magnetism [23-26]. 

For example, the (1 + 3)-splitting of the self-dual repr esentation Zap — Fa/3 + 'i*Fa0 
of a 2-form Fap identifies the complex (spatial) vecto af^)! 

(1) Za = ZafjU^ — Ea+ iBa , ZqU" — , 

where Ea and Ba denote the electric and magnetic parts of the 2-form with respect to 
m", respectively. The relationship between Fap and Za basically establishes an isomor- 
phism between the Lie algebra of the Lorentz group and the Lie algebra o(3, C) of the 
complex orthogonal group in three dimensions. Completing u to an orthonormal frame 
Ea {a — 0,1,2,3) with Eq — u", in terms of which the complex vector Za only has 
spatial components Za {a,b — 1, 2, 3), then Lorentz transformations of this frame result 
in a complex rotation of this complex 3- vector. 

Similarly, the self-dual part of the Weyl tensor (C -t- i*C)ap-yS is associated with the 
complex (spatial) tensor [10] 

(2) - Qa/3 = (C + i *C)a0'ySu''u^ = Eaf) + iBap , 

where now the spatial tensors Eap and Bap denote the electric and magnetic parts of 
the Weyl tensor associated with u", respectively. The spatial tensor Qap is symmetric 
and trace- free, i.e., satisfies the conditions 

(3) (5"q = , Qap = Qpa , Qapu'^ = 

and its nonzero components in a frame adapted to u" are a complex symmetric tracefree 
3x3 matrix {Qat) (a, ^ ~ li2,3) when expressed in an orthonormal frame adapted to 
u". The mixed form of this spatial tensor with component matrix Q = {Q°'b) can be 
thought of as a linear transformation on the representation space of complex 3- vectors, 
and as such its matrix can be classified by its eigenvalues and eigenvectors. 
The eigenvalue problem for this 3-dimensional matrix 

(4) Q\r^ = Xr\ 

is equivalent [10] to the corresponding problem for the Weyl tensor 

(5) ^C'^^sX^^ = AX"'' , 

where X°'f^ are the eigenbivectors of the Weyl tensor associated with the eigenvalue A, 
represented as complex vectors in the self-dual representation. The algebraic type of 

(}) A spatial quantity with respect to is such that any contraction with u" gives 0. 
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the matrix Q provides an invariant characterization of the gravitational field at a fixed 
spacetime point. This Petrov classification based on the eigenvalues Ai,A2,A3 of Q, 
their multiplicities, and the number of linearly independent eigenvectors leads to certain 
canonical forms for the matrix Q called the normal forms which are listed in [10]. At 
most two eigenvalues are independent due to the tracefree condition on the matrix 

(6) TrQ = Ai+A2 + A3 = 0. 

This classification is insensitive to the scale of the eigenvalues, and so only depends on 
their scalefree part. 

The most general Petrov type I corresponds to Q being diagonalizable. A canonical 
or normal frame Ea for this type is one in which the matrix is diagonal 

Ai 

(7) Q/ = I A2 

As^ 

Type D occurs as the subcase Ai = A2 and permutations. The remaining canonical forms 
of Q can be found in [10]. 

The Petrov classification distinguishes between algebraically general spacetimes (type 
I) and algebraically special ones (types D, II, N, III and O) depending on the degeneracy 
of eigenvalues and eigenvectors associated with the above spectral problem. One way 
of determining whether a given spacetime is general or special involves evaluating the 
invariants 

(8) / = ItVQ^ = liXi + XI + XI), J= ^TVQ3 = i(A? + A^ + Af) . 

For Petrov types I and then D and II in the top two levels of the Petrov hierarchy, both 

invariants / and J are nonzero, while both vanish for types O, N and III in the lowest 
level of the three level pyramid. Any algebraically special spacetime satisfies 

(9) = 27 

so except for most degenerate types O, N and III in the Petrov hierarchy, one can intro- 
duce the so called "speciality index" which is a well defined function of the eigenvalues 
Ai , A2 , A3 by taking the dimensionless ratio (which is therefore independent of the overall 
scale of the eigenvalues) 

having the value 1 for the algebraically special types D and II and obviously invariant 
under any permutation of the eigenvalues. 

This quantity is related to the similar invariant quantity M introduced earlier by 
Mcintosh and Arianrhodtf [27-29] by the relation 



(11) 



S = {l + M/27)-'^ . 
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They used this quantity to simphfy the discussion the geometry of the set of principal 
nuU directions of the Weyl tensor described by Penrose and Rindler [30] , and apphed this 
insight to the Kasner and Curzon metrics. 

The tracefree condition Ai + A2 + A3 = does not constrain the overaU scale of the 
three eigenvalues, but it does lead to only one independent ratio among them. Choosing 
/i — A1/A2, then it follows from the constraint that A3/A2 = — (1 + /i) and therefore 
^3/^1 = —(1 + so the six possible ratios ^ab — ^a/\ [a 7^ b) are given by the 

following formulas in terms of one of them n — 

(Mi2, Ai32, Ai3i) f ^ j = (m, ~(1 + Ai), -(1 + , 



, N _ / A2 A2 Ai 

lAi21, Ai23, Ail3j — T-, T-, T- 



Ai A3 A 

(12) (Ai,A2,A3) = A2(m, !,-(! + m))- 

Thus any scale invariant (i.e., dimensionless) function of the eigenvalues can only depend 
on jjL, which must be the case for the speciality index. The invariance of the speciality 
index under the action of the symmetric group 5*3 will reflect these six possibilities for a 
choice of the single independent ratio, and the speciality index will therefore be invariant 
upon replacing fi by any of the five remaining expressions fiab given in terms of /i. 

Using the trace- free property of the matrix Q, one can re-express any eigenvalue in 
terms of the remaining two, so the speciality index is only a function of two independent 
eigenvalues, which we choose to be Ai and A2 

(13) S-'' ^?^^(^^+^^)^ 



4 (A2 + A2 + AiA2)3 



Since the speciality index is a dimensionless homogeneous function of these two variables, 
it only depends on their complex ratio, leading to a convenient parametrization of this 
function 



2 112 



^''^ ''^"^^T (l + ,. + ^^)3 =- [(^+l)2+3l3 - Me 



4J 



S{u) is a complex analytic function of /i (clearly symmetric under the reflection /i + i ^ 

— (/i + ■^)) whose values along the real axis are conflned to the closed interval [0,1], 

1 

which vanishes in the limit /i — )■ 00, and whose poles are the roots /i = — 2±*2 '^^ 
the denominator jj? + ^ + 1, at which values the condition |Ai| = IA2I = IA3I holds. 
Values of /x on the real axis correspond to the vanishing of the magnetic part of the Weyl 
tensor in a canonical frame, the case of a purely electric Weyl tensor. The poles are the 
intersection points of the unit circles |/Lt-|-l| = l = |^|. Indeed these two circles and the 
line through their two intersection points (a line which can be thought of as the circles of 
infinite radius with centers on the real axis at positive and negative infinity) divides up 
the complex plane into six equivalent regions under the action of the symmetric group. 
The imaginary part of S{u) is zero exactly on these two circles and on the line through 
their intersection points plus along the real axis. 

In Figs. [1] and [2] the real and imaginary parts of the speciality index are plotted as 
functions of the real and imaginary parts of /i. In the real part plot, the speciality 
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Fig. 1. - The real part of the speciaUty index function is plotted as a function of the 

real and imaginary parts of The speciality index values corresponding to vacuum Kasner 
spacetimes (lying above the real ax:is) are indicated by a thick solid curve. The Petrov type D 
cases (iS = 1) are denoted by bullets on this curve above /^ = —2, —1/2, 1. 



index values corresponding to purely electric Weyl tensors lie above the real axis, shown 
superimposed as a thick curve, and only half the plot is shown since it has even symmetry 
about the real axis; it also has even symmetry about the line 3?e(/x) = ^1/2. In the 
imaginary part plot, again only half the plot is shown since it has odd symmetry about 
the real axis; it also has odd symmetry about the line 9?e(/i) = —1/2. Below it will be 
seen that the real speciality index over the real part of fi describes exactly the family of 
Kasner spacetime Weyl tensors. 

To get a better feeling for the symmetries of these graphs, one can introduce polar 




Fig. 2. - The imaginary part of the speciality index function S{fj,) is plotted as a function of the 
real and imaginary parts of /i. 
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Fig. 3. - Contours for the values —1,0, 1 of the real (left) and imaginary (right) parts of the 
speciality index function 5 (pi) are plotted using compactified polar coordinates in the fi complex 
plane centered on the point (—1/2,0). The imaginary part shows clearly the two circles |/i + 
1/2| = 1 = l/il and the vertical line through their intersection points, and which are also shown 
in gray in the real part plot, which divide the complex plane up into six equivalent regions under 
the symmetric group. The imaginary part is on these curves and along the real axis. iS(/i) 
has additional reflection symmetries about the center. The dashed curve in the real part plot is 
where the real part is 1, and its intersection with the gray curves 9(5) — give the speciality 
points where S = 1, namely, three points along the real axis. 



coordinates centered on the point (—1/2, 0) and bring infinity in to the unit circle centered 
at that point by compactifying the radial coordinate 

u + 1/2 = tan(r7r/2) exp{i6) . 

Now implicitly plotting the contours where the real and imaginary parts of iS(/i) equal 
-1,0,1 yields Fig.H 

The original eigenvalue permutation symmetries of this function correspond to the 
following changes in u associated with the three basic transpositions 

Ti2(Ai, A2, A3) = (A2, Ai, A3) : 

T23(Ai,A2,A3) = (Ai,A3,A2) : ^l ^ + fr^y^ , 

(15) Ti3(Ai,A2,A3) = (A3,A2,Ai) : /i^-(l + M), 

so S{^) is unchanged under each of these transformations. The fixed points of these 
three transformations are respectively T12: ^ = 1,-1, T23: = 0,-2 and T13: fi — 
— 1/2, = 0, while the two circles |/i + 1| = 1 = and the line /i + 1 = connecting 
their intersection points (the poles) are invariant submanifolds under their action. The 
two cyclic permutations are compositions of the transpositions, explicitly 

^'23i(Ai, A2, A3) = (A2, A3, Ai) : /I ^ -(1 + fi)-^ , 

(16) P3i2(Ai, A2, A3) - (A3, Ai, A2) : M ^ -(1 + fi-') , 
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which both have the same fixed points /j, = — ^ ±i^, which are the poles of <S(/Lt). It 
should be noted that each of these transformations of the complex plane are combinations 
of inversions, reflections and translations. 

The special types D and II have two eigenvalues which are the same, corresponding 
to one of the three eigenvalue ratios being 1, related to each other by the permutation 
symmetries 



(17) 



A2/A1 = 1 
A3/A1 = 1 
A3/A2 = 1 



M=l, (Ai,A2,A3)cx(l,l,-2), 
M = -2, (Ai,A2,A3)a(-2,l,l), 
/x = -l/2, (Ai, A2, A3) cx (1,-2,1), 



so 1 = S{1) = S{—2) = <S(— 1/2). The zeros of S are obviously /z = 0, — 1 and in the 

limit corresponding to a third zero at infinity, again all related to each other 

by the permutation symmetries. These correspond to another privileged triplet of values 
which occur when two of the eigenvalues are of opposite sign 



(18) 



A2/A1 = -1 
A3/A1 = -1 
A3/A2 = -1 



M = -l, (Ai, A2, A3) (X (1,-1,0), 
M = 0, (Ai, A2, A3) cx (0,1,-1), 
M-^=0, (Ai, A2, A3) cx (1,0,-1). 



3. " The Kasner family 

The Kasner metric is 



(19) 



ds^ = -dt^ + t^P'dx^ + t^P-dy^ + t^P'dz^ 



where the so-called Kasner indices (constants) satisfy 



(20) 



Pi +P2 +P3 = 1 = Pl +P2+P3 



SO that they form a one-parameter family of triplets which can assume values in the 
interval [— ^,1]. The space of distinct 3-gcomctries corresponds to the quotient of this 
family of spacetimes by the discrete group of permutations of the spatial coordinates 
and can be represented by any one of the six intervals of ordered values of the triplets: 
Pa < Pb < Pc, (a, b, c) = cr(l, 2, 3), cr e 6*3. The Kasner family admits two distinct special 
subcases for which two of the Pa indices are equal: 1) Pa = Pb = 0, Pc = ^ for which 
the spacetime is flat and therefore of Petrov type O and 2) pa — —1/3, Pb = Pc = 2/3 
corresponding to the Kasner locally rotationally symmetric type D solution, with 5=1. 
The matrix of mixed spatial extrinsic curvature tensor components {K°-b) = —l/2{g°-'^gcb) 

h 1/2 

-^''a 111(566 )■ is diagonal 



(21) 
(22) 
(23) 



-t"Miag(j3i,p2,P3) , 

-t~^{pi +P2+Pz) = 

diag(pi,p2,P3), 



from which it is clear that the Kasner indices represent the scale free part of this tensor's 
eigenvalues. The Kasner index constraints come from the scale invariant Hamiltonian 
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constraint in vacuum with zero spatial curvature 

(24) iK\K\ - K^.K\)/iK%r = + fa + g + Ps)^ ^ ^ 

{pi+P2 + P3r 

which forces the corresponding paths in the minisuperspace of inner products to be null 
geodesies of the deWitt metric, and from the trace of the spatial Einstein equations, 
which is an evolution equation for the square root of the spatial metric determinant 

^1/2 _ ^Pl+P2+P3 = f 

(25) = g-^/^{k\y = g-'/^{g'/^y= (pi + p2 +Ps)(pi + P2 + Ps - 1)^"' • 

Of course these same two independent Einstein equations are also satisfied by pa = 0, 
which is just flat spacetime, not so interesting. 

Expanding the square and cube of the constraint pi + P2 + P3 = 1 and simplifying 
using both constraints leads respectively to two other identities 

1 1 1 ^^ 
P1P2 + P2P3 + P1P3 = 0, (^ 1 1 =0) 

Pi P2 P3 

(26) pI+pI+P3 = '^ + ^PlP2P3 ■ 

The first (which is the scale invariant combination of the two Kasner constraints) shows 
that the three products of the Kasner indices (equivalcntly their reciprocals) parametrize 
the three diagonal entries of a diagonal tracefree matrix modulo the overall scale. Similar 
using the first of these two identities to replace piP2 in their product and then using the 
linear constraint to eliminate pi + p2 leads to 

(27) P1P2P3 = -^3(1 - P3) or P1P2 = -P3(l - Ps) • 

There are a number of parametrizations of this family of Kasner indices by a single 
independent parameter. For example, the obvious procedure is to solve the basic pair 
of contraints in terms of any one of the indices, say pi. Eliminating p^ using the linear 
constraint leads to a quadratic equation in p2 with coefficients depending on pi , whose 
obvious solution (backsubstituting to express pa in terms of pi as well) [7, 8] is 

m = ^ ([1 - Pi] ± [(1 - + 3pi)]^/^) , P3 = 1{[1- Pi] T [(1 - Pi)(l + 3pi)] V2) . 

On the other hand the spatially flat Kasner spacetime plays an important role as a 
limiting case in the dynamics of spatially homogeneous spacetimes with more complicated 
symmetry, which gives an extra context for choosing the parametrization conveniently. 
The very different trigonometric parametrization quoted by Stcphani et al [10] turns out 
to be related to Misner's approach to this dynamics, while the one introduced by Lifshitz 
and Khalatnikov describes extremely efficiently a particular situation in this dynamical 
arena. 

The Lifshitz-Khalatnikov parametrization [13-16] represents the Kasner indices in 
terms of a single independent (real) parameter u = U32 = P3/P2 as follows 

-u 1 + u u{l + u) 

(2^) = \ 2 ' P'^ = T, \ 2 ' P^ = \ 2 • 
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-1/3 



2/3 







Fig. 4. - The Kasnor indices {pi,p2,pa) arc plotted versus u, together with the speciality index S. 
The obvious reflection symmetry across u = —1/2 shows 3 pairs of corresponding intervals in this 
parametrization, but hides the equivalence of all 6 intervals corresponding to the permutations 
of Pa < Pti < Pc- Pi and p2 and the speciality index S all have the limiting value at ±oo while 
P3 goes to 1. 



For a given ordered set of three distinct Kasner indices, there arc six permutations that 
correspond merely to permuting the spatial coordinates, dividing the real u axis up into 
six equivalent intervals (see Fig. 1) whose endpoints are the set {— oo, —2, —1,-1/2, 0, 1, oo} 
Confining u to the interval u E [l. oo) leads to the ordering of the Kasner exponents: 
Pi < P2 < P3- In general Pa < Pb < Pc for Ucb = Pc/Pb G [1, oo). 

The Kasner index parametrization appears mysterious until one thinks about the two 
constraints on the indices. Both constraints fix the overall scale of the indices so that 
only their ratios can vary, so it is natural to pick one of those ratios as a new parameter, 
which can be done in six different ways. Picking u = U32 = P3/P2, for example, it is 
trivial to express all six possible ratios in terms of it. Eliminating = up2 first, and 
then using the linear constraint to eliminate pi = 1 — P2 — up2 = 1 — (1 + u)p2, the 
quadratic constraint reduces to a quadratic equation in u 



with roots P2 = (obvious) and p2 = u/{l + u + u'^), from which pi and ps are then 
obtained by backsubstitution. One can define six such parameters Uab = Pa/Pb, for which 
Pc l£ Pb l£ Pa when (c.b.a) is a cyclic permutation of (1,2,3), all of which are needed 
for application to the various pairs of axes in BLK-Mixmaster dynamics. This is quite 
unremarkable, but on the other hand this simple class of parametrizations turns out to be 
quite remarkable in the BLK-Mixmaster dynamics, as will be discussed in detail below. 

By straightforward evaluation [4,7], one finds that the complex matrix (Q"b) is already 
diagonal, real and purely electric (since the magnetic part of the Weyl tensor is zero), 
with the value 



(30) 



[1-(1 + «)P2]'+P2 + («P2)' = 1 



(31) 



(g%) = -{E\) = -K\ + K\K\ 

= -t"^diag(pi - p\,p2 - pI,P3 - pI) 
= t~'^diag{p2P3,P3Pi,PiP2) , 
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SO Kasner is in general of Petrov type I and its eigenvalues are 

(32) (Ai,A2,A3) = t^'^(p2P3,P3Pl,PlP2) , 



so that 



,oQ^ ^2 P3 A3 P3 /-, , X Ai P2 , _i, 

33 — = — =u, — = — = -l + tt, — = — = -{l + u ), 

A3 P2 Ai pi A2 Pl 

hence /i23 = u is purely real for these spacetimes. The permutations of the eigenvalues 
are locked to those of the Kasner indices, which were originally discussed by Misner 
[17] using logarithmic diagonal metric variables adapted to the pure trace and tracefree 
decomposition of diagonal matrices 

(ln(5i?),ln(522^'),ln(.9^3/') = W\ [3^ , P') 

(34) = 1, 1) + /3+(l, 1, -2) + /?- V3(l, -1, 0) . 

The Kasner indices are then functions of the ratios of the corresponding cotangent space 
(momentum) variables (po,P+,P-)- 

To visualize the permutation symmetry in a symmetric way, it is useful to parametrize 
the indices using an angular polar coordinate in the diagonal tracefree matrix cotangent 
space following Stephani et al [10] 

1 1 1\ 2 / ^ , 27r, , 47r, 

Pi - ^,P2 - g,P3 " 3 ) = 3 ( cos^',cos(^' + y),cos(5' + y) 

(35) = -i[cos^'(-2,l,l) + sin«'\/3(0,l,-l)], 

o 

where here a permutation of Misner's tracefree diagonal matrix parametrization (|34p 
occurs, here adapted to the first axis rather than the third axis as usual. Fig. [5] shows 
the Kasner indices plotted as the signed radial variable versus the polar angle 4* measured 
counterclockwise from the half ray along the direction rotated by 7r/3 counterclockwise 
from the positive horizontal axis. The triplets marking each of the six sectors give the 
ordering of the Kasner indices for that interval. 



4. BKL-Mixmaster dynamics 

The Kasner indices themselves or the parameter u or the angular parameter ^' may 
be thought of as parametrizing the unit circle in the p+-p- plane of the two anisotropy 
momenta which coordinatize the cotangent plane to the tracefree metric matrices, them- 
selves parametrized by the Misner [17] variables and /3~, as explained in detail in [20]. 
Letting (j)^,p_) be the restriction of these variables to the unit circle, one has 

{p+,p^) = (cos(* + 7r/3), sin(^' + 7r/3)) 
(36) = -{u^ + U+ iy\[u + Ijlf - 3/4, \/3[u + 1/2]) . 



Since the DeWitt metric on the diagonal metric matrices is conformally Minkowskian 
(proportional to — (i/3o^+(rf/3'*')^ + (rf/3~)^, one can identify the unit circle in the anisotropy 




Fig. 5. - The Kasner indices {pi,P2,P3), identified respectively by short dashed, long dashed 
and solid cardiods, are plotted in the "Misner plane" of the p± momentum variables along the 
radial direction versus the translated polar angle 3" measured counterclockwise from the half 
ray at an angle 7r/3 in the counterclockwise direction from the positive horizontal axis, where 
the negative Kasner indices lead to the internal loops. This establishes a correspondence of the 
Kasner indices and the parameters u and ^ with the unit circle of directions in this p+-p- plane, 
whose axes are adapted to the third direction E3. The three circles have radii 1/3, 2/3, 1, while 
each of the six sectors is labeled by a triplet (o, 6, c) indicating the ordering Pa < P6 < Pc in that 
sector. 



momentum space with the unit circle in the anisotropy velocity space, representing the 
space of directions in the anisotropy plane. This "Kasner circle" represents the 
asymptotic Kasner states in a scattering off a curvature potential in spatially homoge- 
neous dynamics idealized as a collision with a moving potential wall in this anisotropy 
plane. 

As explained in detail elsewhere [20] , a collision with each of the three potential walls 
arranged with equilateral triangle symmetry is idealized as an exact Bianchi type II 
solution for the diagonal metric components. In terms of the u = U32 Lifshitz-Khalatnikov 
parametrization of the Kasner circle, with the future time asymptotic value Mqo € (1, 00) 
and past time asymptotic value u_oo € (—00, —1), a collision corresponds to the system 
point moving from wall 2 towards wall 3 while colliding with wall 1 . The simple geometric 
reflection which takes place in the rest frame of the moving wall is boosted when described 
in the nonmoving frame. 

The pictorial representation of the Kasner map in terms of a collision with a curvature 
potential wall in this anisotropy plane is complicated by the fact that the dynamics is 
reversed in time (as one approaches an initial cosmological singularity, the cosmic time 
t decreases) while the potential wall is moving, so the analogy is more like pinball than 
billiards since the moving wall acts like a flipper in transferring momentum to the system 
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Fig. 6. - The Kasner circle divided into six zones by its scattering properties with curvature 

walls W1,W2,W3 parallel to the sides of the equilateral triangle. The L-K parameters M31 
and 'U32 describe the successive collisions of the system point with Wl, W2 as it moves into 
and out of the corner C3. The momentum scattering can be imagined as taking place at the 
origin. For a case in which the system point is moving from wall W2 towards wall W3 while 
colliding with wall Wl as time increases, if the final state is in zone (123): u € (1, 00), the initial 
state is in zone (231)U(213): u e (—00, —1), with a collision map which is simply u — » —u. In 
the increasing time direction, wall Wl is moving towards the corner CI with half the speed 
of the system point, so even if the system point is moving at an angle up to 30 degrees from 
the direction of the wall Wl (dashed axis) in zone (213), a collision still occurs. The Kasner 
dynamics goes in the reverse time direction. 



point. Fig. 6 indicates the orientations of the three equilateral triangle curvature walls of 
BLK-Mixmaster dynamics relative to the Kasner circle (not the walls themselves, which 
move in the [3^ plane, see Fig. 7. 

Each of these curvature walls moves with a speed which is half the speed of the system 
point itself, allowing a relative motion of the system point of up to 30 degrees away from 
a direction parallel to a straight line wall and still collide with the wall as the walls move 
inward as the cosmic time increases. For example, in Fig. 6 the incoming unit vector 
with values in regions (2,3,1) and (2,1,3) ranging from direct normal incidence (u = —1) 
to moving away at an angle 30 degrees from the normal (u = — oo) is kicked by the wall 
into the region (1,2,3). 

Of course it is not enough to look at the diagram in velocity/momentum phase space. 
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Fig. 7. - The anisotropy plane interior to the three curvature walls W1,W2,W3 forming 
an equilateral triangle is divided into six zones by the perpendicular bisectors of the sides, 
conveniently labeled by the corresponding six Kasner circle intervals. In turn these pair up into 
three corner zones CI, C2, C3. The L-K parameters U31 and U32 of the previous figure describe 
the successive collisions of the system point with Wl, W2 as it moves into and out of the corner 
C3. 



One must couple this with the actual position in the f3 plane relative to the actual 
triangular walls. See Fig 7. The three bisector line segments of the equilateral triangle 
angles extending from the vertex to their common meeting point at the center divide up 
the triangle into three isoceles triangles, each of which has one of the three walls Wl, 
W2, W3 as a base, while extending those bisectors through the center to the opposite 
side further divides these into six right triangles, one for each sector of the Kasner circle. 
These in turn can be matched up in pairs corresponding to the three corners CI, C2, 
C3. Fig. 6 then describes the directions in the tangent space to a system point located 
in one of these six triangles. 

Suppose the system point is in corner C3 with a final direction in the Kasner interval 
(1,2,3), namely M32 = u G (1,00), so that it is escaping towards a collision with wall 
W3 after a collision with wall Wl. Its initial incoming direction was in the Kasner 
interval (2,1,3), coming from a collision with wall W2, whose final direction corresponds 
to = —u, U31 = u — 1. If M — 1 > 0, then the initial state for this collision will have 
U31 = —{u — 1), U32 = u — 2. Working backwards in time, we are describing a Kasner 
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era consisting of a series of incoming-outcoming Kasner states (epochs) which continue 
until u — n < 0, and the system point escapes corner C3 moving on to a Kasner era 
bouncing between the two walls of another corner. To graphically record this dynamics 
in the /3^-plane, one rescales the wall position variables in order to fix them in place, 
giving an equilateral triangle billiard representation of the more complicated moving wall 
bouncing of the system point, thus removing the scale in anisotropy space. See the figures 
in Chapter 13 of Ryan and Shepley [18], taken from Ref. [19]. 

The transitions in the scalefree part of the Weyl tensor eigenvalues during each col- 
lision with a curvature wall are remarkably simple in terms of the L-K parameters. For 
example, if it > 2 then in the direction towards the initial singularity a collision with 
wall Wl in corner (73 (the standard way in which the Ksner map is usually stated for 
ordered Kasner indices) leads to 

(Ai,A2,A3) = A2(u, 1,-1-u) A2(-u, 1,-1-l-u) 
(37) = A'i(l,u- 1,-u). 

The first eigenvalue ratio is just reversed in sign, and the final ratio compensates to keep 
the sum zero, with respect to the same eigenvector basis. Transposing the basis as well 
leads to the unit subtraction of the old ratio to form the new ratio appropriate for the 
new basis. 

5. — Concluding remarks 

While not a very fundmental revelation, finally the mysterious Lifshitz-Khalatnikov 
Kasner index parametrization is understood in a larger context than the narrow one in 
which it was introduced without explanation in order to do its job. It is satisfying to 
pull on a bit of mathematics and discover that it is hiding a lot more than meets the 
eye initially. Many of us were attracted to relativity in part for its mathematical beauty, 
and these little gems are always a joy to discover. 
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